home *** CD-ROM | disk | FTP | other *** search
/ SGI Freeware 2002 November / SGI Freeware 2002 November - Disc 3.iso / dist / fw_plotutils.idb / usr / freeware / share / ode / dynamo.ode.z / dynamo.ode
Text File  |  2002-01-08  |  682b  |  41 lines
  1. # This is a disk dynamo simulation of the earth's magnetic field.
  2. # It attempts to simulate the dramatic switch in the polarity
  3. # of the field that takes place every eon or so.
  4. #
  5. # Source: Michael Steele, B.A. thesis, Reed College, 1981 (Physics)
  6.  
  7. # You may run this example by doing:
  8. #
  9. #    ode < dynamo.ode | graph -T X -C
  10. #
  11. # or alternatively, to get a real-time plot,
  12. #
  13. #    ode < dynamo.ode | graph -T X -C -x 0 10 -y -10 10
  14. #
  15.  
  16. # The equations are: 
  17. #
  18. # w' = Q -zy - Vw
  19. # y' = S( Az - y)
  20. # z' = wy - z
  21. # S = 5
  22. # V = 1
  23. # Q = 14.625
  24. # A = 1
  25.  
  26. A = 1
  27. V = 1
  28. Q = 14.625
  29. S = 5
  30.  
  31. w' = Q - z * y - V * w
  32. y' = S * ( A * z - y)
  33. z' = w * y - z
  34.  
  35. w = 1
  36. y = 1
  37. z = 1
  38.  
  39. print t, z
  40. step 0, 10, .01
  41.